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Abstract. We present [C i] 3 Pi- 3 Po spectra at four spiral arm positions and the nuclei of the nearby galaxies M83 
and M51 obtained at the JCMT. The spiral arm positions lie at galacto-centric distances of between 2kpc and 
6kpc. This data is complemented with maps of CO 1-0, 2—1, and 3-2, and ISO/LWS far-infrared data of [Cll] 
(158 fim), [O i] (63 /im), and [N n] (122 fim) allowing for the investigation of a complete set of all major gas cooling 
lines. From the intensity of the [Nil] line, we estimate that between 15% and 30% of the observed [C n] emission 
originate from the dense ionized phase of the ISM. The analysis indicates that emission from the diffuse ionized 
ly-j ■ medium is negligible. In combination with the FIR dust continuum, we find gas heating efficiencies below ~ 0.21% 

in the nuclei, and between 0.25 and 0.36% at the outer positions. Comparison with models of photon-dominated 
regions (PDRs) of Kaufman et al. (1999) with the standard ratios [O l](63)/[G h]pdh and QO i](63) + [C ii]pdr) 
vs. TIR, the total infrared intensity, yields two solutions. The physically most plausible solution exhibits slightly 
lower densities and higher FUV fields than found when using a full set of line ratios, [C ii]pdr/[C i](1-0), [Ci](1- 
0)/CO(3-2), CO(3-2)/CO(l-0), [C ll]/CO(3-2), and, [Oi](63)/[Cii] PD r. The best fits to the latter ratios yield 
densities of 10 4 cm -3 and FUV fields of ~ Go =20-30 times the average interstellar field without much variation. 
At the outer positions, the observed total infrared intensities are in perfect agreement with the derived best 
fitting FUV intensities. The ratio of the two intensities lies at 4-5 at the nuclei, indicating the presence of other 
mechanisms heating the dust. [Cl] area filling factors lie below 2% at all positions, consistent with low volume 
filling factors of the emitting gas. The fit of the model to the line ratios improves significantly if we assume that 
[Cl] stems from a larger region than CO 2-1. Improved modelling would need to address the filling factors of the 
various submm and FIR tracers, taking into consideration the presence of density gradients of the emitting gas 
. by including cloud mass and size distributions within the beam. 

Key words. Galaxies: ISM, Galaxies: structure, Galaxies: individual: M83, M51, ISM: Structure, Infrared: Galaxies, 
Submillimeter 



1. Introduction 

Neutral atomic carbon is thought to form predominantly 
in surface layers of molecular clouds where C n recom- 
bines and CO is dissociated due to the far-UV photons 
governing the chemical reactions. FUV photons (6eV< 
hv < 13.6 eV) are primarily responsible for the heating of 
the surface regions via photoelectric effect on dust grains 
while at larger depths cosmic-ray induced heating will 
dominate. These regions are referred to as photo dissoci- 
ation regions o r, more generally, as photon dominated re- 
gions (PDRs) ijTielens fc Hollenbachlll985t IStoerzer et all 
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Il996t iKaufman et alJll999|) . PDR models take into ac- 
count the relevant physical processes, solve simultaneously 
for the chemistry (using an extensive chemical network) 
and the thermal balance, as a function of cloud depth. It 
is found that the r atio of [C III / fC i] is an accurate tracer 
of the FUV field (jCerm fc Phillips! l2000h . parametrized 
by Go in unit s of the Habing-field 1.6 10~ 3 erg s _1 cm -2 
Another important parameter governing 
the depth at whi ch C I forms is the ratio of density over 
FUV field n/G l|Tielens fc Hollenbachlll985|l . This ratio 
also determines the efficiency of converting FUV photons 
to gas heating, i.e. the photoelectric heating efficiency e 
l|Ra.kes fc TielenJlTflfl^. 
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While the Milky Way survey o f FIR lines conducted 
with COBE/FIRAS l|Fixsen et al.ll999j) showed that [C n] 
is the dominant cooling line, it also showed the importance 
of the two finestructure lines of [C i] . Both lines are ubiqui- 
tous and the two lines together amount to 75% of the total 
cooling of all rotational CO lines in the inner galaxy. This 
picture has also emerged from extragalactic observations 
of the [Ci] 1-0 line. These show again that the cooling 
due to Ci and C O are of the same order of magnitude 
for most ga laxies ((Bavet et al.ll2004t llsrael fc Baasll2003l 
I2002L 120011) . Ci is found to be a good tracer of molecu- 
lar ga s, possibly more reliable than CO l|Gerin fe Phillips 
l2000l) . 

Several coordinated mapping studies of nearby galax- 
ies have been started during the past years. The BIMA 
SONG survey l|Regan et al.ll200l"h has aimed at obtaining 
the 12 CO emission of 1-0 and 2-1 rotational lines at high 
spatial resolutions. There exist velocity-integrated [C 11] 
observations of lar ge samples of g alaxies with the KAO 
llStacev et a,l.lll99lh and with TSO llMalhotra, et a,lJl200ll) 
herea fter MKHOl, and (|Leech et alJll999t INegishi et alJ 
120011) at ~ 1' resolution. Th e SINGS Spitzer Legacy 
Project ijKennicutt et alJ l2003[) has started imaging 75 
galaxies in the infrared, including M51. In the coming 
years, both SOFIA and the Herschel Space Observatory 
are expected to provide velocity-resolved [C 11] data at res- 
olutions of ~ 10", complementary to many current single 
dish observations of CO and [Ci]. 

In external galaxies where a large number of clouds 
or even GMCs fill the beam it is difficult to seperate 
the different contributions and judge their importance. 
A long standing problem is that a substantial fraction 
of the [C ll] emission may originate from the diffuse ion- 
ized and neutral medium. Comparison with [Nil] helps 
to estimate the fraction originating from PDRs, but usu- 
ally with large uncertainties due to the varying chemi- 
cal and excit ation conditions in d ifferent galactic envi- 
ronments fe.g. lContursi et al"l l2002'). The present study is 
part of the preparatory work for future airborne and space 
missions like SOFIA and Herschel. In addition, it may 
serve as template for studies of e.g. [C 1] and CO in high-z 
galaxies which have recently become possible JWeiss et al.l 
SlWalter et al.l2004tlPetv et a,lJ2004lNeri et a,l.l2003t 
Weiss et al.ll2003|) . 

Most extragalactic observations of atomic carbon have 
so far concentrated on the bright galactic nuclei or en- 
hanced emission of edge-on galaxies. Here, we compare ob- 
servations of the two nuclei of M83 and M51 with pointed 
observations at spiral arm positions which show enhanced 
star forming activity. The galacto-centric distances of the 
selected outer positions lie between 1.8 and 5.8 kpc. We 
combine observations of atomic carbon with low and mid- 
J CO and 13 CO data, as well as FIR [Cn], [Oi](63^m), 
and [Nil] (122 (im) data from the ISO data-base and thus 
include the brightest gas cooling lines of the far-infrared 
and submillimeter regime. 



Table 1. Basic properties of M83 and M51. D25 x (fes is 
the optical diameter from the RC3 catalogue. Lpir is the 
surface integrated far-infrared luminosity between 42.5 and 
122.5 /im using the listed distance. L60/L100 is the IRAS color 
ratio of the surface integrated flux densities at 60 /jm and 
100 fim, corrected for extinction. References: (1) RC3 cata- 
logue ^^^uco^^^^^l Ej|yl), (2) iTalbot et alJ dl979l): 
iTullvl Jl98sl). (3) iRice et al.l Jl988h . (4l lSandage fc Tammannl 
^Irie VaucoukuTsl fl97^ 





M83 


M51 


RA(2000) 


13:37:00.5 


13:29:52.7 


DEC(2000) 


-29:51:55.3 


47:11:43 


Type 


SAB(s)c (1) 


SA(s)bc pec (1) 


Distance [Mpc] 


3.7 (5) 


9.6 < 4 > 


10" correspond to 


179 pc 


465 pc 


Heliocentric velocity [kms -1 ] 


516« 


463 (1) 


Position Angle [deg] 


45 < 2 > 


170 


Inclination [deg] 


24 (2) 


20 


-D25 x d 2 5 ['] 


12.9 x 11.5 (1) 


11.2 x 6.9 (1) 


Lpir [10 9 L ] 


7J (3) 


14(3) 


-F60/-F100 


0.43 (3) 


0.44 (3) 


Leo [Jy] 


286 (3) 


85< s > 



1.1. M83 

M83 (NGC 5236) is the most nearby CO-rich grand-design 
spiral galaxy, seen almost face-on (Table^J. It has a pro- 
nounced bar, with two well-defined spiral arms connected 
to the starb urst nucleus. In this paper, w e adopt a distance 
of 3.7 Mpc (jde Vaucouleurs et al.lll99l]l though recent ob- 
servations of Cepheids indica te a slightly larger distance 
of 4.5 Mpc l|Thim et al.ll2003h . 

Lo w - J CO maps we r e obtain ed bv ICrosthwaite et alJ 
ll2002l): iLundgren et alJ i|2004bllah : iDumke et all ll200lh : 
ISakamoto et alJ l|2004tl . [C i] observations of th e cen- 
ter were conduc t ed b y llsrael fc Baas! l|2002j) and 
IPetitpas fc Wilson! l|l998h . Pointed KAO observations re- 
port strong FIR fine-structure lines tow ards the nucleus 
with a rapid fall-off towards the arms ^Crawford et alJ 
Il985h . 

Here, we present new [Ci] data of the center and 
two spiral arm positions on the north-eastern arm and 
south-western bar-spiral transition zone. The emission of 
[Cn], [Nil] (122), and [O in] (88) observed with ISO/LWS 
(Brauhcr 2005, priv. com.) is strongly enhanced in these 
interface regions indicating greatly enhanced star forma- 
tion rate s. ISO/LWS emissio n from the center was ana- 
lyzed bv INegishi et al . (2001). The north-e a stern arm 
was previously st udied bv lLord fc Kennevl l)l99l|) and 
IRand et alJ {l999) who presented OVRO interferometric 
12 CO 1-0 maps. The eastern position at (89", 38") pre- 
sented h ere correspo n ds to the bright feature #6 com- 
piled bv IRand et all l|l999l Table 4). About 15" to the 
east of the CO arm newly formed stars form the op- 
tical arm and an HI ridge. At #6, the CO and dust 
arms coincide while they are offset further to the south. 
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Position (—80", —72") studied here corresponds to a CO 
1-0 peak in the south-western bar-spiral transition 
zone which exhib its a massive GMC co mplex and lumi- 
nous Hn regions ijKennev fc Lordlll99lj) . Note that only 
less than 5% of the s ingle dish flux is recovered by the 
interferometric maps ijRand et al.l fl999). Thus, relatively 
smoothly distributed diffuse molecular gas is completely 
missed. 



1.2. M51 

The nearby grand-design spiral galaxy M51 (NGC 5194) 
seen almost face on (Tabled) is interacting with its small 
companion NGC 5195, which lies 4.5 ' to the north. M51 
is a Seyfert 2 galaxy l|Ho et alJll997h. The central A GN 
is surrounded by a ~ 100 pc disk ( Kohn o et al ] ITflfl3) of 
dense and warm gas ^Matsushita et alJll998IK 

A large amount of observational data are available 
for this object, inc l uding an extended KA O map of [C n] 
l)Nikola et al.ll200lh . lOa.rnett et all l|2004l) used ISO/LWS 
data of the M51 H n region CCM 10, to study the ionized 
gas, finding that abundances are roughly solar. 

Several single-dish studies mapped the low-lying ro- 
tational 12 CO and 13 CO transitions. CO 1-0 and 2- 

lo et al] l|l993bO and 



1 was mapped by 
iNakai et all Jl994li : 

design spiral, CO 
arms. Maps of CO 



Garcia-Buri 



iTosaki et all ( 2002). In this grand- 
is tightly confined to the spiral 
3-2 and 4-3 were obtained at th e 
; IWielehinski eTaTI lll999h : 
who show that warm molecular gas 



HHT hvlNieteu et a,lJ (1 
iDumke et all |2£ 

is extended in M51 at galacto-centric distances of at least 
upto 100", resp. 5kpc. Single-dish observations of neutral 
carbon were so far obt a ined only in the cent er region by 
iGerin fc Phillips! l|2000h : llsrael fc Baasl l|2002|b and Israel, 
Tilanus, Baas, 2005, in prep.). 

Ap e rture synthesis maps we re obt ained bvlAalto et alJ 
l|l999|) : ISakamoto et~ai1 l|l999l) and iRegan et alJ l|200lj) 
CO 1-0 at reso lutions of 4" - 6". Recently, 



iMatsushita et alJ l(2004h mapped the inner region in CO 
3-2. 

Here, we selected two positions at the spiral arms 
lying in the northeastern and the southwestern zones, 
i.e. at 72", 84" and -84", -84", of e nhanced [C III emis - 
sion tracing enhanced star formation l|Nikola et alJl200ll) . 
ISO/LWS data are available fo r these positions, and for 
the center iNegishi et alJ ()200l|) an d Brauher (200 5 . priy . 
com.). The Hil region studied bv iGarnett et alJ l)2004|) 
using ISO/LWS, CC M 10, lies at about +148",+45" 
jCarranza et alJll969|b 1.5' to the north-west of 72",84". 



2. Data sets 

We present here observations of [Ci], CO, and 13 CO spec- 
tra at four spiral arm positions and the centers of M51 
and M83 (Tabled . We combine these with ISO/LWS FIR 



Table 2. Overview of the [Ci], CO, and 13 CO data of 
M83 and M51. 9b is the telescope half power beamwidth. 
All spectra are scaled to main beam brightness tempera- 
tures, T m b = T^/t) m b. References: 1: new data for this pa- 
per; 2: ICrosthwait e"et alJ <l2002ft ; 3: Bayet et al. (priv. com.); 
4: IWielebinski et alJ lll999t):lNieten et al.1 Jl99^ : lDumke et alJ 
fcOOlT) : 5: iGarcia-Burillo et alJ Jl993bl) . The last column in- 
dicates whether single pointings or maps are available. The 
mapped CO data was gauss-smoothed to the 80" ISO/LWS 
resolution. 



Line 


e h 

(") 


Telescope 


»7mb 


Refs. 




M83 at (0,0), ( 


89,38), (-80, -72) 








[Ci] 1-0 


10 


JCMT 15m 


0.52 


1 


point. 


CO 1-0 


21 


IRAM 30m 


0.78 


1 


point. 


CO 1-0 


55 


NRAO 12m 


0.88 


2 


map 


CO 2-1 


10 


IRAM 30m 


0.57 


1 


point. 


CO 2-1 


28 


NRAO 12m 


0.56 


2 


map 


CO 3-2 


25 


CSO 10m 


0.75 


3 


map 


13 CO 1-0 


22 


IRAM 30m 


0.79 


1 


point. 


13 CO 2-1 


10 


IRAM 30m 


0.60 


1 


point. 



M51 at (0,0), (72,84), (-84,-84) 








[Ci] 1-0 


If) 


JCMT 15m 


0.52 


1 


point 


CO 1-0 


21 


IRAM 30m 


0.65 


5 


map 


CO 2-1 


10 


IRAM 30m 


0.46 


5 


map 


CO 3-2 


22 


HHT 10m 


0.50 


4 


map 


13 CO 1-0 


22 


IRAM 30m 


0.79 


1 


point. 


13 CO 2-1 


10 


IRAM 30m 


0.60 


1 


point. 



spectral line data at all six positions together with the FIR 
continuum derived from HIRES/IRAS 60 /im and 100 /im 
maps. 

2.1. [Ci] data taken at the JCMT 

We have observed the fine structure transition of atomic 
carbon ([Ci]) at 492 GHz (609 ^m, 3 Pi- 3 P ; hereafter 
1—0) in M51 (3 positions) and M83 (3 positions) using the 
JCMT 15m telescope. We used receiver RxW with a sin- 
gle mixer and the DAS autocorrelator. Observations were 
carried out during 35 hours in May and June 2003. We 
used the double-beamswitch observing mode with a wob- 
bler throw of ±3' in the direction of the major axis, i.e. 
in the direction of the largest velocity gradient. Pointing 
was checked using SCUBA after an initial alignment with 
RxW and Jupiter at the start of each shift. It was found to 
be accurate to within 2 — 3". The atmospheric zenith opac- 
ity at 225 GHz varied slowly between 0.1 and 0.05, corre- 
sponding to a r of 2 and 1 at 492 GHz. After merging the 
DAS autocorrelator spectra using the SPECX software, 
further data analysis was done using the CLASS/GILDAS 
package of IRAM. 

2.2. CO data taken at the IRAM 30m MRT 

We have observed the 12 CO and 13 CO 1-0 and 2-1 rota- 
tional transitions at all six positions in M51 and M83 using 
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Fig. 1. Spectra of M83 at the central and two spiral arm 
positions (Table[2J). Offsets are given in arcseconds relative 
to the (0,0) position (Tabled . Velocities are relative to 
LSR. All spectra are on the T m b scale. The 12 CO data are 
at a common resolution of 80". The [Cl] and 13 CO data 
are at their original resolutions listed in TableE] 

the IRAM 30m telescope. These observations were carried 
out in double beam switch mode with a wobbler throw 
of ±4' using the filterbank of 1 MHz resolution for the 
3 mm band and the 4 MHz filterbank for the 1 mm band. 
Observations were carried out during 15 hours on July, 
23rd and 26th, and on September, 10th, 2004. Pointing 
and focus were checked and corrected every ~ 2 hours. 
The pointing accuracy was better than 4". The amount 
of precipitable water vapour varied slowly between 10 and 
~ 20 mm. Telescope parameters are listed in Tabled 



;ions in the spiral arms of M83 and M51 




300 400 500 600 300 400 500 600 400 500 600 
Velocity (km/s) Velocity (km/s) Velocity (km/s) 



Fig. 2. Spectra of M51 at the central and two spiral arm 
positions (cf. Fig.^). 

Table 3. FIR continuum and line intensities in units of 
10~ 6 ergs -1 cm' 2 sr " 1 . The absolute calibration error is 
assumed to be 15% l|Grv et al.ll200l| . 



Positions 


HIRES 


ISO/LWS 






(Aa, AS) 


FIR 


[Cn] 


[Nn] 


[Oi] 


(",") 




158/mi 


122/im 


63/j.m 


M83: 










(0,0) 


32 10 3 


81.8 


14.1 


86.7 


(-80,-72) 


7.2 10 3 


38.6 


5.3 


27.3 


(89,38) 


6.6 10 3 


33.4 


5.9 


26.6 


M51: 










(0,0) 


17 10 3 


44.1 


12.3 


32.2 


(72,84) 


4.6 10 3 


16.7 


< 2.3 


13 


(-84,-84) 


3.5 10 3 


15.4 


2.1 


13.4 



2.3. FIR line fluxes from ISO/LWS 

The central area of M83 covering ~ 5' x 4' was mapped 
on a fully-sampled grid of 61 positions with ISO/LWS. 
M51 was observed at 13 positions, mainly along a cut 
through the center and the two prominent [C n] lobes 
in the north-east and south-west seen in the KAO map 
hv lNikola etall (|200lh - The ISO/LWS line emission data 
was uniformly processed by Brauher (2005, priv. com.) 
to derive line fluxes in Wm~ 2 . To convert to intensities, 
we use a LWS beam size of 80" (SIlws = l-210~ 7 sr), 
the mean value published in the latest LWS Handbook , 
and extended source corrections factors llOrv et alJl200l . 
Resulting ISO intensities of [C n] (158 /xm, [Nil] (122 /urn), 
and [Oi] (63 fim) at the positions observed in [Cl] are 



listed in TablcOH Since the [O i] (146 /im) line was detected 
only at the center positions (S. Lord, priv. comm.), we did 
not include it in the present analysis. 

2.4. FIR continuum maps taken with IRAS 

To derive the far-infrared continuum at all selected po- 
sitions, we obtained high-resolution (HIRES) 60 /im and 
100 /im IRAS maps from the IPAC data center 1 . Enhanced 
resolution images were created after 200 i terations using 
the maximum correlation method (MCM lAumann et alJ 

1 A description of IRAS HiRes reduction is available at 
http : //irsa. ipac . caltech. edu/IRASdocs/hires_over .html 
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Il990l) . These data were smoothed to an effective com- 
mon circular beam of 80" and then combined to create 
m aps of the fa r -infra red flux. The FfR flux is defined as 
in lHelou et all Jl988l) : FIR = f .26 KT 14 [2.58 7(60 /im) + 
1(100 ^m)] where FIR is in Wm~ 2 and I is in Janskys. 
FIR is a good estimate of the flu x contained between 42.5 
and 122.5 ^m l|Helou et al.ll988|) . TableEJlists the FIR flux 
at the selected positions. To derive the total i nfrared flux 
TIR, we follow the procedure introduced by iDale et ail 
(2001) who derived an analytical expression for the ratio 
of total infrared flux TIR to the observed FIR flux, i.e. the 
bolometric correction, as a function of the 60/im/lOO/im 
flux density ratio from modelling the infrared SEDs of 
69 normal galaxies 2 For M83 and for M51, the total in- 
frared flux is a factor 2.3 larger than FIR given the global 
60 ^m/100 /im ratio of 0.43 (Table©). 



3. Spectra and line ratios 

The main aim of this work is to use the com- 
bined FIR ISO/LWS, HIRES/IRAS, and the millime- 
ter/submillimeter line data coherently and to investigate, 
to what degree these give a consistent fit within the frame- 
work of a simple model scenario such as PDR excitation. 
We therefore smoothed the 12 CO maps (Table© to the 
ISO/LWS angular resolution of 80" using Gaussian ker- 
nels. Spectra of [Ci], CO, and 13 CO are displayed in 
Figures lTl2*l and integrated intensities in Tabled Ratios 
with [Ci] and 13 CO for which no maps exist, were cor- 
rected for beam filling (Tabled . 

3.1. M83 

The [Ci] lines are widest at the center position with 
130 kms" 1 FWHM and drop to about SOkms" 1 at the 
two outer positions. In addition, peak line temperatures 
drop strongly, leading to a pronounced drop of [C i] inte- 
grated intensities and area integrated [C i] luminosities by 
factors of 10 to 18 at galacto-centric distances of less than 
2kpc (Table©. 

The CO 3-2 transition traces warm and dense gas since 
its upper level energy corresponds to 2.8 J( J+ 1) = 33.6 K 
and the critical density needed to thermalize this line is 
4 10 3 J 3 ~ 1 10 5 cm~ 3 , only weakly dependent on the ki- 
netic temperature. Trapping typically reduces the critical 
densities by up to an order of magnitude, depending on 
the optical depth of the lines. The CO 3-2/1-0 line ratio 
of integrated intensities (Table© thus is a sensitive tracer 
of local densities for densities of less than ~ 10 5 cm -3 . 
Here, both line maps were smoothed to 80" resolution, 
beam filling factors thus cancel out to first order. The es- 
timated calibration error is 21%. While the 2-1/1-0 ratio 
is ~ 1 at all positions indicating that the J = 2 state is 
thermalized, the 3-2/1-0 ratio drops slightly from 0.60 in 



2 Note that TIR corresponds to the bolometric FIR dust 
continuum emission I(FIR) used in the PDR models of 



iKaufman et al.1 iT99' 



the center to 0.46 at the spiral arm positions. These ra- 
tios indicate that densities are lower than 1 10 5 needed for 
thermalization of the J — 3 level. 

The 12 CO 2-1 vs. 13 CO 2-1 line ratio, as well as the 
corresponding 1-0 ratio, trace the total column densities 
of the 13 CO line. The ratios remain constant at ~ 10 for 
2-1 and ~ 9 for the 1-0 transition. 

The above ratios found at the center position agree 
within the quoted errors with the ratios presented in 
llsrael fc Baasl (|200lh . 

In M83, we find a significant drop of the CI vs. 13 CO 

2- 1 ratio from about 4 in the center to ~ 1.4 at the two 
outer positions. 

3.2. M51 

Again, the [Ci] line is widest at the center, ~ 100 kms -1 
FWHM, and drops to ~ 30kms _1 at the two outer posi- 
tions. The 13 CO lines show similar line widths at all three 
positions. Peak temperatures of [Ci] hardly drop between 
the center and the two outer positions. Resulting inte- 
grated intensities and luminosities drop by a factor of ~ 4 
only. 

The observed CO 2-1/1-0 ratios lie between 0.6 and 
0.8 for all three positions and do not peak at the center. 
This indicates that not even the 2-1 line is thermalized 
in M51. However, the CO 3-2/2-1 ratios do not drop as 
would be expected but equal the CO 2-1/1-0 ratio or even 
exceed them, while staying below 0.8. This is difficult to 
explain with a single component model as we will show 
below. The center CO 3-2/1-0 ratio of 0.55 is in agreement 
with ratios previously found with single-dish telescopes 
which range between 0.5 — 0.8 at beam sizes of ~ 14" 
to 24" ^Matsushita, et allll999l iMauersberger et alJll999t 
IWielebinski et al.lll999() . Interferometric observations at 
~ 4" resolution trac ing the dense nucl e ar ga s show a high 

3- 2/1-0 ratio of 1.9 lMatsushita et all l|2004h . 

High central column densities are indicated in this 
work by the rather low 12 CO/ 13 CO 2-1 ratio of 4.6. In 
contrast, the outer postions show ratios between 6 and 
14. 

We find CI vs. 13 CO 2-1 ratios of ~ 1 at the center 
and at (-84, -84) while (72,84) exhibits a high ratio of 
3.2. 



3.3. Comparison 

The gradients of [C i] luminosities with galacto-centric dis- 
tances are strikingly different in M51 and M83. M83 is 
much more centrally peaked, [C i] luminosities drop by a 
factor of 18 at only 1.8 kpc distance in M83. In contrast, 
luminosities in M51 drop by only a factor of 4 at galacto- 
centric distances which are more than a factor of 3 larger, 
i.e. at 5.8 kpc. However, the central [Ci] luminosities of 
M51 and M83 agree within 40%. 

The line widths observed at the outer positions of M51 
and M83 are typical for the disks of these two galax- 
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ies llHanda et all Il99(j iGarcia-Burillo et all Il992|) . See 
Table 2 of iGarcia-Burillo et alJ l|l993a) for a compilation 
of CO line widths found in these and several other galax- 
ies. 

The CI vs. 13 CO 2 -1 ratio in t he Mil ky Way is of - 
ten found to b e 1 (e .g. lKeenelll995|) while llsraell l|2004(l : 
llsrael fe Baasl l|2002fl find a strong variation of this ratio 
for 15 galactic nuclei. The [C i] line is stronger than the 
13 CO 2-1 line for all but three galaxy centers. The highest 
ratios are about 5. Here, we find a variation between 1.3 
and 4. 

For galaxy centers, llsrael fe Baasl ()2002|) found that 
this ratio is well correlated with the [C i] luminosity cov- 
ering a range of 160Kkms _1 kpc 2 in the active nucleus 
of NGC3079 down to ~ 1 Kkms -1 kpc 2 in the quiescent 
center of Maffei2. Here, we increase the range down to 
0.11 Kkms -1 kpc 2 (TableE|) at the same resolution of 10". 
In contrast to the galaxy centers, the spiral arm positions 
observed here do not show a systematic correlation be- 
tween the [Ci]/ 13 CO line ratio and [Ci] luminosity. 



4. Physical Conditions 

4.1. Simple homogeneous models 

As shown above, the observed CO line ratios cannot be 
explained with a simple LTE analysis. As a first step in 
order to estimate the kinetic temperatures and local den- 
sities of the CO emitting gas, we present in this section the 
results of slightly more realistic escape probability radia- 
tive transfer calculations of ho mogeneous spherical clumps 
dStutzki fc Winnewisserll985l). We assumed a 12 CO/ 1 3 CO 
abundance ratio of 40 l|Mauersberger fc Henkell 19931) . As 
input, we use the four ratios of integrated intensities 
(Table©: 12 CO 3-2/2-1, 12 CO 2-1/1-0, 12 CO/ 13 CO 1- 
0, and 12 CO/ 13 CO 2-1. We calculated model intensi- 
ties of the three transitions for column densities 10 14 < 
A^(CO)/Av/(cm- 2 /(Kkms- 1 )< 10 22 , local H 2 densities 
10 < rtioc/cm~ 3 < 40, and kinetic temperatures 10 < 
2kin/K< 40. The modelled ratios were compared with the 
observed ratios taking into account the observational er- 
ror of 21% to derive the \ 2 ■ The best fitting iV(CO)/Av, 
n, and T^in together with the corresponding minimum re- 
duced x 2 are listed in TableHO 

Below, we first describe the results of the one- 
component fits for the positions observed in M83 and M51. 
Much more complete physical models of the emitting re- 
gions are presented in the next sectior JOl 

M83: At 80" resolution, the 12 CO line ratios do not vary 
significantly between the center and the two bright spiral 
arm positions observed here. A 12 CO 3-2/1-0 ratio of 0.5 
indicates an excitation temperature of ~ 10 K assuming 
optically thick thermalized emission and simply using the 



detection equation. The escape probability analysis leads 
to a similar result (Table©. At all three positions, the 
ratios are well modelled by a one-component model with 
a rather low kinetic temperature of only 12 — 15 K and 
a density of rti oc = 310 3 cm~ 3 . This result does not ex- 
clude the existence of a warmer and denser gas phase as 
would be traced by higher CO transitions or the 63 iim 
[Oil li ne as discussed below. For the center. llsrael fc Baasl 
(2001) have in fact deduced a warmer phase by including 
observations of the CO 4-3 line in their radiative transfer 
analysis. 

The 12 CO and 13 CO line ratios found in M83 are char- 
acteristic for dynamica lly active or starburst region s in the 
classification scheme of lPapadopoulos et al . (2004). In this 
scheme, extreme starbursts would show similar 12 CO 2- 
1/1-0 and 3-2/1-0 ratios but much higher 12 CO/ 13 CO 
ratios. 

M51: In the escape probability analysis of CO and 13 CO 
ratios in M51, we discarded solutions leading to temper- 
atures below 12.5 K and densities below 10 3 cm -3 as un- 
physical. Since the CO J =3-2 is strong at all positions, 
densities and temperatures must be higher. 

The CO 3-2/2-1 ratios are ~ 0.8 at all positions 
in M51, significantly higher than in M83. This indicates 
higher densities than found in M83. Indeed, the escape 
probability analysis finds densities between 3 10 3 and 
310 4 cuT 3 . The 12 CO/ 13 CO ratios agree with this solu- 
tion for low temperatures of ~ 12 K. However, the ob- 
served CO 2-1/1-0 ratios are too low to agree with this 
solution. Neither the temperatures nor the densities are 
well constrained, and minimum chi squared values are 
high. This shows the short coming of a one-component 
model even when trying to model only the three lowest 
rotational CO tra nsitions. And it is in fact in agreement 
with the finding of lGarcia-Burillo et alJ l|l993bj) who used 
the lowest two transitions of 12 CO and 13 CO and couldn't 
find a set of Tkin and ni oc fitting simultaneously the line 
ratios for the arms and for the central position. 
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Table 4. Integrated intensities in Kkms- 1 (all columns but col. 1, 2 and 4) at the resolutions used for the spectra in 
Figs. 11121 Column 2 gives galacto-centric distances i? ga i in kpc. The calibration error is estimated to be ~ 15%. (The 
observ ed [ C i] intensities at the cent er positions agree within 20% with previous data presented in lGerin fe Phillips! 
( 2000) and lPetitpas fc Wilson! (|l998).) Column 4 lists [C i] luminosities in Kkms 1 kpc 2 in brackets, i.e. the intensities 
integrated over the 10" JCMT beam. 
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Table 5. Line ratios of integrated intensities in Kkms 1 . The la errors are ~ 21%. The CO 3-2/1-0 ratio was added 
to allow an easier comparison with the literature. To compare all data at a common resolution of 80", we scaled the 
[Ci] and 13 CO data, for which maps are not available, using the beam filling factors <i>^ ^ 10 and <i>^ ^ 21 . The factor 
$ 8 ° /w for [Ci] and 13 CO 2-1 data was derived from the ratio of 12 CO 2-1 integrated intensities at 80" and at 10" 
resolution. Likeweise, the factor <&^ 21 for 13 CO 1-0 data was derived from CO 1-0. 
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Table 6. Results of the escape probability analysis using the four observed line ratios 12 CO 3-2/2-1, 12 CO 2-1/1-0, 
12 CO/ 13 CO 1-0, and 2-1. JV(CO) is the total CO column density p er 80" beam. M is th e total mass integrated over 
the beam and assuming an abundance ratio of [CO]/[H 2 ]=8.5 10~ 5 l|Frerking et alJll982|) . n av is the average density 
assuming that emission stems from a sphere with beam diameter. 
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4.2.1. [Cn] emission from the ionized and neutral 
medium 




Fig. 3. Comparison of the observed intensity ratios 
[O i] (63 /um)/[C ii W and([Oil(63 um) + fC n] PDR )/TIR c 
with PDR models ijKaufman et aljfl999(l at the six posi- 
tions in M83 and M51. 



4.2. PDR Analysis 

To further constrain the physical conditions at the ob- 
served positions in M83 and M51, we compare the ob- 
served line intensity ratios with the resul ts of the model 
for Ph o ton Dominated Region s (PDRs) bv lKaufman et all 
l|l999|) : lTielens fc Hollenbachl l|l985|) . The physical struc- 
ture is represented by a semi-infinite slab of constant 
density, which is illuminated by FUV photons from one 
side. The model takes into account the major heating 
and cooling processes and incorporates a detailed chemi- 
cal network. Comparing the observed intensities with the 
steady-state solutions of the model, allows for the deter- 
mination of the gas density of H nuclei, n, and t he FUV 
flux, G n, where Go is measured in units of the me 
l)l968|) value for the average solar neighborhood FUV flux, 
1.6 x 1 0~ 3 ergs cm" 2 s" 1 . As has been pointed out in de- 
tail bv iKaufman et aTl (fl999) in their Section 3.5.1 and 
several other authors, the application of these models to 
extragalactic observations is not straightforward since in- 
dividual molecular clouds are not resolved in single-dish 
observations and several phases of the ISM are therefore 
observationally coexistent within each beam. The addi- 
tional many degrees of freedom in the parameter space 
for more complex models, however, are ill constrained by 
the few observed, beam-averaged line ratios. Hence, sim- 
plistic models, e.g. with only a single source component, 
are used to at least derive average properties of the com- 
plex sources. Nevertheless, it is possible to obtain some 
insight into the spatial structure and the local excitation 
conditions, as we will show. 



Carbon has a lower ionization potential (11.26 eV) than 
hydrogen, so that [C n] emission arises not only from pho- 
ton dominated regions, but also from the ionized phases 
of the ISM and from the diffuse neutral medium traced by 
Hi. 

Analyzing the Milky Way FIR line d at a ob - 
tained with FIRAS /C OBE l|Fixsen et alJ Il999j) . 
IPetuchowski fc BennettJ (|l993l) argue that slightly 
more than half the [Cil] emission in the Milky Way 
arises from PDRs, the remainder from the extended low 
density warm ionized medium or diffuse ionized medium 
(DIM), and an insignificant portion from the ordinary 
cold neutral medium (CNM). 

[Cil] emission from the ionized medium. Here, we use the 
observed [Nn](122^m) and [Cil] lines, to estimate the 
fraction of [C n] originating from the ionized medium. 
However, a thorough analysis would need more FIR 
emission line data from the ionized medium, in par- 
ticular the [Nii](205/im) line, to discriminate the rel- 
ative importa nce of the differen t phases of the ion- 
ized medium ((Bennett et al.lll99^ . Extragalactic obser- 
vations of the [N ill (205/n m) are however very rare to date 
l|Petuchowski et alJli99^) . 

The components of the ionized phase of the ISM which 
contribute to the [C n] emission are dense HII regions 
(n e > 100 cm" 3 ) and the d iffuse ionized medium (DIM) 
(n e - 2cm" 3 ) l|Heileslll994 . 

The fraction of [Cil] stemming from Hn regions de- 
pends st rongly on th e elect ron density of the ionized 
medium. ICarral et alJ l|l994h showed that upto 30% of 
[C n] stems from H II regions when electron densities 
exceed 100cm -3 . For dense HII regions (n c ^> n cr ), 
model calcula tions suggest [Cii]/[Nii] = 0.28(C/N) dcnsc 
teubirJll985|) where (C/N) dcnso is the abundance ratio. 
The Galactic abundance ratio found in dense Hn regions 
is 3.8 l|Rubin et alJll98SHl993h . We thus expect to find an 
intensity ratio of 



[Cll]/[Nll] ionjH n= 1.1. 



(1) 



Observations of radio free-free emission suggest that 
the DIM has typical volume densities of n c ~ 2 cm" 3 
Mezeerl ll97S ) and temperatures of T c = 4000 K 
Mueller et al. Il987|) . Densities are thus considerably 
lower than the critical densities of 310 cm -3 for the 
[N III (122) line and 50 cm" 3 for the [Cn](158) line l(Genze< 
1991). In this limit, the total power emitted in the col- 
lisionally excited fine-structure lines is simply the pho- 
ton energy times the upward collision rate, as given 
by iHeiled (^994). The resulting ratio of intensities is: 
[Cii]/[Nii](122) = 3.05(C/N) diff . Galactic absorption 
line measurements of diffuse gas give (C/N)dim = 1-87 
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Table 7. Physical para meters at the observed positions M83 and M83, derived from fitting the observed intensity 
ratios to PDR models of iKaufman et al.1 l)l999|) . Column (2) gives the galacto-centric distance. Columns (3) and (4) 
list the fitted local densities and FUV fields together with the corresponding surface temperature (col. (5)) and the 
minimum chi squared (col. (6)). The filling factor (f> vv — G , o bs/G reflects the ratio of TIR C to the best fitting FUV 
field G - 0a 1 is the 

area filling factor of the [C i] emitting regions, i.e. the ratio of observed [C i] intensity vs. that of 
the best fit model, corrected for beam and velocity filling (see text). 
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M51: 


















(0,0) 


0.00 


4.0 


1.25 


66. 


3.2 


1.88 


4.31 


0.013 


(72, 84) 


5.40 


4.0 


1.25 


66. 


5.2 


1.32 


1.17 


0.016 


(-84, -84) 


5.77 


4.0 


1.25 


66. 


14.1 


1.19 


0.88 


0.009 



( Sofi a et al Jl997tlMever et alll997ft . which implies an ex- 
pected intensity ratio of 



([Cii]/[Nii](122)) ionjD iM = 5.7. 



(2) 



The procee ding section follows bas ically the arguments 
of MKH01 and lContursi et all l|2002l) . 

The metallicities, parametrized by the Oxygen abun- 
dances, have be en found to be only s lightly supersolar in 
M83 and M51. IZaritskv et ail l|l994|) find (O/H) abun- 
dances of 1.5 10~ 3 and 1.9 10~ 3 respectively, at 3kpc 
galacto-centric distance, from visual spectra of Hn re- 
gions, which is about a factor 3 higher th an the solar 
metallicity of 0.46 10~ 3 l|Asplund et al.ll2004| ). Abundance 
gradi ents with radius are found to be shallow in M83 and 
M51 l)Zaritskv et al Jl994h . More recent observations with 
IS O/LWS and modelling of the CCM10 H II region of M51 
bv lGarnett et all |2004) indicate instead that the (O/H) 
abundances are about a factor of 2 less, i.e. roughly solar. 
In addition, l|Garnett et al.lll999(l showed that the (C/N) 
abundance ratio, which is of interest for the [Cii]/[Nii] 
ratio discussed here, is independent of metallicity in both 
normal and irregular galaxies. We therefore use Galactic 
abundances to estimate th e inten s ity ra tios in M83 and 
M51. For the Milky Way, iHeilej l)l994l) have estimated 
that [Nil] originates predominantly from the DIM, con- 
tributing ~ 70%. This was derived from the observed 
[Nll](122)/[Nii](205) ratio using photo ionization models. 
However, the [Nil] (205) line has not yet been observed in 
M83 and M51. We therefore use Eqs.EJand^ assuming 
that [Nil] stems solely from Hn regions, or alternatively, 
solely from the DIM. 

Next, we can then derive the fraction of [C n] emission 
originating from PDRs: 



[Cll] 



PDR 



[Cn 



obs 



(ft!L x [ n "« 122 >- 



(3) 



The ratios of observed [C n] versus [N n] intensities 
vary between 3.6 and 7.4 at the six positions observed 



in M83 and M51 (Tabled ■ IcTarnett et all l|2004h found a 
higher ratio of 11.8 at the CCM 10 position in M51. 

The ratios found in M83 and M51 lie at the low end of 
the ratios MKH01 found in the sample of 60 unresolved 
normal galaxies studied. They show a mean ra t io of 8 and 
a scatter between 4.3 and 29. IContursi et al.l (j2002) ob- 
served ratios of more than 7.7 in N GC6946 and ratios o f 
greater than 4 and 10 in NGC1313. Ifligdon et all l|2003l) 
found ratios between > 2.6 and 20 in M33. 

If the [Nil] emission originates only from the diffuse 
ionized medium, then the major fraction of [C n] arises 
from this phase, and only a small fraction from PDRs 
(Table®. The observed [Cii]/[Nii] is <5.8 at 3 positions, 
including the two nuclei, which would indicate that no 
[C n] emission at all arises from PDRs. This is however un- 
realistic, since the emission of the [Oi](63/im) line, stem- 
ming from warm, dense PDRs, is strong compared to the 
[Cn] lines, especially in the nuclei (Table®. For this rea- 
son, we discard this solution. 

Assuming on the other hand, that the fraction of [C n] 
from the ionized medium and all [N n] emission stem only 
from the dense Hn regions, then a fraction of only 15% 
to 30% of [Cn] originates from this phase (Eq.®, while 
70% to 85% of the observed [C n] emission then stems 
from PDRs. We prefer this solution and use it in the PDR 
analysis discussed belo w. In an ISO/LWS study of star- 
forming regions in M33,|mid on et alJ l|2003|) have recently 
used FIR lines of the ionized medium, i.e. [O in](88/im), 
[O m](52/mi) and others to estimate the electron densities 
and other parameters of the emitting gas, estimating that 
between 7% and 47% of [C n] stems from H n regions, for 
their sample of positions. They conclude that the DIM is 
not needed to explain the observations. 



[Cn] emission from the diffuse neutral medium. The pre- 
dicted [C n] emission from the atomic gas has in general, 
for many galactic nuclei, been found to be far too weak 
to account for the observed [Cn] emission l(Stacev et alJ 
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Table 8. FIR line ratios after correcting the observed 
[C n] intensities for emission from the ionized medium. We 
study two cases: The emission of [Cn] from the ionized 
medium originates entirely a. from the diffuse medium, 
and, b., that the emission arises entirely from dense Hn 
regions. Also shown is the intensity ratio e of the two ma- 
jor cooling lines of PDRs ([Oi](63)+[Cii]pdr) vs. the ob- 
served TIR continuum. 
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0.81 
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0.21 


-80,-72 


7.3 


0.85 


0.83 


0.36 


89,38 


5.7 


0.81 


0.99 


0.35 


M51: 










0,0 


3.6 


0.69 


1.06 


0.16 


72,84 


7.3 


0.85 


0.91 


0.25 


-84, -84 


7.4 


0.85 


1.02 


0.33 



because the density is not high enough to appre- 
ciably excite the [Cn] emission at the meas ured Hi col- 
umn densities. This view was confirmed by l|Carral et alJ 
Il994|) who conducted a detailed study of FIR cooling 
lines of NGC253 and NGC3256. Both in M51 and in 
M83, no large-scale correlation between Hi emission and 
that of [C n] is seen, indicating aga in that [Cil] does 
not trace the diffuse neutral medium ijNikola et al.ll200ll 
ICrawford efall|l985l). 

iNikola et al.l ll200ll) used H I col umn densities 
( Tilanus & Alien! Il99lt iRots et all Il99(t to derive the 
contribution to the [Cil] emission in M51, assuming the 
same range of temperatures, densities, and ionization 
fractions for the warm and cold neutral medium (WNM, 
CNM) as have been found for the Milky Way. They find 
that the contribution of the WNM is negligible for most 
of the M51 disk except the northwest, which was not 
studied here. The contribution of the CNM is estimated 
to be less than 10%— 20% in all regions but the northwest. 

We have thus not corrected the [C 11] emission for a 
possible contribution from the diffuse neutral medium. 

4.2.2. The infrared continuum and the FUV field 

The stellar FUV photons heat the molecular gas and 
dust which subsequently cools via the FIR dust contin- 
uum and, with a fraction of less than ~ 1% (|Stacev et alJ 
ll99lh .MKH01 via [Cil], [Oi](63), and other cooling lines. 



To the extent that filling factors are 1 and other heating 
mechanisms like cosmic ray heating can be neglected, the 
observed TIR continuum intensity should equal the mod- 
elled FUV field. This is also expected if a constant fraction 
of FUV photons escape without impinging on cloud sur- 
faces. 

The PDR model of iKaufman et al ] lll999h assumes 
a semi-infinite slab illuminated from one side only. For 
the extragalactic observations described here, we however 
have several PDRs within one beam and the clouds are il- 
luminated from all sides. Hence, the optically thin total IR 
intensity stems from the near and far sides of clouds. Here, 
this i s taken into account by dividing the observed TIR 
by 2 (Kaufm aTet alll999h : TIR C = TIR/2 = 2.3 x FIR/2 
(cf. Sec l2.4(l . While this correction holds exactly only for 
finite plane parallel slabs illuminated from both sides, it is 
a good first approximation. The corrected TIR can then 
be used to derive the corresponding FUV intensity via 
G ,obs = TIR r 47r/ (2 x 1.6 10~ 3 ) er g s cm~ 2 s" 1 . Following 
the arguments of IKaufman et all l(l999h . the additional 
factor 2 takes into account equal heating of the grains 
by photons outside the FUV band, i.e. by photons of 
hv < 6cV. We find a variation by one order of magni- 
tude, 15 < G 0:O bs < 144 (Tabled Fig.0}. 

4.2.3. First estimates of FUV field and density 

The two intensity ratios [Oi](63)/[Cii]pdr and 
([Oi](63)+[Cii] PD r)/TIR c , of the two major PDR 
cooling lines and the continuum, have been used ex- 
tensively to derive the density and FUV field of the 
emitting regions (e.g. MKH01). Since [Oi](63) and [Cil] 
are the dominant coolants, the latter ratio is a good 
measure of the phot oelectric heating efficiency e (e.g. 
IKaufman et al.l ^999). The former ratio measures the 
relative importance of [Cil] vs [Oi](63) cooling. For high 
FUV fiel ds and high densities, the ratio becomes larger 
than one l|Kaufman et al.l f l9 99). 

In M83 and M51, the intensity ratio [O i](63)/[Cii] PD r 
varies only slightly between 0.8 and 1.3 (TableJHb)- The 
heating efficiency varies between ~ 0.25 and 0.36% at the 
outskirt positions while it drops to below 0.21% in the 
centers. 

The values which we find in M83 and M51 lie within 
the range covered by MKH01, who find heating efficiencies 
ranging between ~ 0.3% and ~ 0.05% for the 60 galaxies 
studied, while the [O i]/[C 11] ratios range between 0.3 and 
~ 10. Though the scatter is large, the heating efficiency 
tends to be high > 0.15% for [O i]/[C 11] ratios of less than 
2. The [Oi](63)/[Cii]pdr ratios found in M83 and M51 
agree roughly with the average value found by MKH01, 
while the heating efficiencies in M83 and M51 span the av- 
erage value of MKH01 upto the highest efficiencies found 
by them. 

MKH01 corrected the observed [Cil] emission by 
roughly 50% when taking into account the contribution 
from the ionized medium, based on the Milky Way re- 
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suits. Here, we corrected by only 15% to 30% (Table|SJ}). 
This uncertainty in how best to correct the [C II] fluxes, 
needs to be considered when comparing the derived heat- 
ing efficiencies and [Oi]/[Cii] ratios. 

The small scatter of the observed two ratios at the 
6 positions in M83 and M51 indicates that the emitting 
gas has similar physical properties. Comparison with the 
results of the Kaufman PDR model shows that two so- 
lutions exist (FigurelHJ. The data can be reproduced ei- 
ther by high FUV fields at low densities or by low FUV 
fields and high densities. The high-Go solution indicates 
2.5 < log(n/cm" 3 ) < 3.2 and 2.4 < logG < 3. As we will 
show, the low-Go solution is less plausible. It indicates 
rather high densities of 4.3 < log(n/cm~ 3 ) < 4.5 and low 
FUV fields of 0.1 < logGo < 0.5. In this case, the observed 
[C il] intensities are more than three orders of magnitude 
larger than the modelled intensities which would indicate 
that many PDR slabs along the lines of sight. Since the op- 
tical depth of the [C ill line in the line c enters is expected 
to be about one (|Kaufman et alj [l999). this scenario is 
discarded. This argument also holds when velocity filling 
is taken into account, since the velocity filling factor is 
< 40 at all position s, as discussed below. We note that 
iHitrdon et all (|2003h in their analysis of ISO/LWS data of 
M33, also discussed the two possible PDR solutions, and, 
following a different line of reasoning, also preferred the 
high-Go solution. 

The n, Go values we find for the high-Go solution, 
agrees with the average values found by MKH01 who also 
exclude the low-Go solution. Their sample of 60 unre- 
solved galaxies covers a slightly larger range of values: 
2 < log(n/cm~ 3 ) < 4.25 and 2.5 < logG < 5. Our values 
also agree with the a verage value found in NGC 6946 by 
IContnrsi et all J2002h. 



the best fit models thus identified. The [Ci](l-0)/CO(3- 
2) and the CO(3-2)/(l-0) ratios are excellent tracers of 
the local densities, almost independent of Go for Go > 10. 
The sensitivity of [C n] on the FUV field is reflected by 
the line ratios which include [C n] . The combination of all 
these ratios therefore allows, in principle, to deduce n and 
Go- 



Densities. Though the quality of the fit varies strongly be- 
tween the six positions, the best fitting FUV field and den- 
sity are almost identical: n = 10 4 cm -3 and 18 < Go < 32 
(TableEJ). Model surface temperatures are ~ 70 K. This is 
slightly lower than the upper energy level of [Cn](158), 
E up /k — 92 K, and much lower than the corresponding 
level of [Oi](63), E up /k = 228 K. Since the critical den- 
sities of the two [O i] lines for collisions with H are also 
high, > 10 5 cm~ 3 , the [Oi] lines are subthermally excited. 
Nevertheless, gas cooling is dominated by [Cn] and [Oi], 
contributions from H2 or Si, covered by ISO/SWS, are 
negligible at these low temperatures. The best fitting den- 
sity agrees within a factor of 3 with the local densities 
derived above via simple radiative transfer analysis, using 
only the 12 CO and 13 CO line ratios (Table©. 



UV filling factor. At the four spiral arm positions, the total 
infrared continuum agrees perfectly with the best fitting 
FUV field (Table[7J), which is an independent confirmation 
of the validity of the PDR analysis using the five line ra- 
tios. The active center regions of M83 and M51 show filling 
factors of ~ 4, indicating that other sources than PDRs 
heat the dust leading to the high observed TIR intensities, 
e.g. massive protostars. 



4.2.4. PDR model fitting of five line ratios 

In order to determine with greater confidence the 
densities and UV fluxes which can explain the in- 
tensity ratios, we have performed a x 2 fitting of 
the observed 5 line intensity ratios [Cii]pdr/[Ci](1- 
0), [Ci](l-0)/CO(3-2), CO (3-2) /CO (1-0), [Cn]/CO(3- 
2), and, [O l](63)/[ C iiIpdr relative to th e predictions of 
the PDR model bv lKaufman et alJ l|l999f) : 

1 / t?obs nmod x o 

2 = 1 \^( th _j V 

X 4-2 <n J 

i—1 

with the observed and modelled ratios R° hs and R^ 1 ^, 
and the error <7j which is assumed to be 20% for all ra- 
tios 3 . In Fig. bleach panel shows the observed line inten- 
sity ratios and the calculated reduced \ 2 i n greyscale for 
each observed position in M83 and M51. Table summa- 
rizes the local density (n) and FUV radiation field (Go) 
at the position of minimum reduced x 2 corresponding to 

3 There are two degrees of freedom since we use four inde- 
pendent line ratios to fit two parameters. 



Quality of the fits. The minimum reduced chi squared of 
the 4 independent ratios lie between 3 and 15 at all po- 
sitions (Tabled). While the two nuclei and the NE-spiral 
arm position of M51 show x 2 of better than 5, the other 
three spiral arm positions show rather poor \ 2 values of 
13-15. 

Inspection of Figure^] shows that, at the latter three 
positions, the [C ii]pdr/[C i] ratio indicates higher FUV 
fields than the best fitting solution. In addition, Figure^] 
also shows that the CO 3-2/1-0 line ratios indicates sys- 
tematically lower densities than the observed [C i]/CO 3-2 
ratios. This holds to varying degrees for all positions and 
for Go > 10. By assuming [Ci] intensities which are a 
factor 2 higher, the quality of the fit is considerably im- 
proved, while the best fitting Go and density stay constant 
at all positions. At (—80,-72) in M83 for example, the % 2 
is improved from 12.9 to 3.4. This indicates that the beam 
filling factors of the [C 1] emission, derived from the CO 2- 
1 data (Table|SJ), are too small, i.e. [Ci] is more extended 
than CO. 

We also conclude that our results are consistent with 
the assumption that only the dense ionized medium con- 
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tributes to the [Cn] emission. There is no need for an 
extended diffuse component. 

Overall, the observed five ratios cannot be well fitted 
with a single plane-parallel PDR model of constant density 
at any of the positions. This is not surprising as the ex- 
citation requirements of the various tracers collected here 
vary widely. Especially, the critical densities vary between 
5 10 2 cm~ 3 for the lower [C i] transition and 5 10 5 cm -3 
for the [O i] 63fim transition. Any density gradients in 
the emitting medium may thus lead to the above discrep- 
ancies with a single PDR model depending also on the 
chemical and temperature structure. In addition, the large 
difference between local densities derived here and beam 
averaged densities of more than three orders of magni- 
tude (TableEJ) shows that the emitting volume must be 
filled with very small but dense structures. From many 
Galactic observations, it is expected that these structure 
show a spectrum of masses, adding to the complexity ig- 
nored here. 



Absolute intensities. When comparing the observed inten- 
sities with the model results, the velocity filling has also to 
be taken into account. The observed [C i] line widths Av b s 
range between 30 and 130 kms -1 FWHM (cf. Sec 1X11 and 
I3.2|l measuring the dispersion of clouds within the beam 
for these extragalactic observations. On the other hand, 
the microturbulent velocity dis persion 5„ of the gas o f one 
PDR model is set to 1.5 kms -1 l(Kaufman et alll99 9). cor- 
responding to a Gaussian FWHM Av mo( j of 3.5 kms -1 as 
is typical for individual Galactic clouds. To calculate the 
[C i] area filling factors (Tabled) > wc divided the observed 
[C i] intensities by the predicted [C i] intensity from the 
best fitting model, corrected for beam (TableEJ) and ve- 
locity filling, viz., 



([ci] obs x *r° 



) Av n 



[CI] 



mod 



Av obs 



[C i] emission fills only a few percent of the 80" beam, 
the derived area filling factors vary between (1 — 2) 10 -2 
(Table[7|), consistent with the low volume filling factors 
described above. 

5. Summary 

We have studied all major submillimeter and far infrared 
cooling lines together with the dust total infrared contin- 
uum at the center positions of the two galaxies M83 and 
M51 and at four spiral arm positions. 

We observed [C i] 1-0 at the six positions at 10" resolu- 
tion. Complementary [Cli], [Oi](63), and [Nil] (122) data 
were obtained from ISO/LWS at 80" resolution. CO maps 
of the lowest three transitions were obtained from the lit- 
erature and smoothed to the ISO/LWS resolution. We also 
obtained pointed 13 CO 1-0 and 2-1 data at all positions. 
In order to allow a comparison of all these data, the [C i] 
and 13 CO data were scaled with beam filling factors de- 
rived from the 12 CO data. For completeness, we also ob- 



tained the total far-infrared continuum intensities from 
HIRES/IRAS 60 ^m and 100 /im data. 

— Integrated intensities peak in the two centers. 
However, M83 is much more centrally peaked than 
M51 as seen in CO and [Ci] as is already seen in the 
spectra fFigs. lTl^ . This is seen even more drastically 
in the drop of [C i] luminosities with galacto-centric 
distance. In M83, luminosities drop by more than one 
order of magnitude over ~ 2kpc, while in M51, they 
drop by only a factor of ~4 over ~ 6 kpc. Obviously, 
this analysis should be refined by more observations of 
[Ci] at different radii. 

— The 12 CO 3-2/1-0 line ratios lie below 0.6 at all posi- 
tions, indicating subthermal excitation of the 3-2 line, 
i.e. densities are less than 10 5 cm -3 . This is confirmed 
by more detailed analysis using escape probability and 
PDR models. The former homogeneous models indi- 
cate densities between 3 10 3 and 3 10 4 cm -3 from CO 
and 13 CO line ratios. 

— We estimated the fraction of observed [C n] emission 
originating from the ionized medium by using the [N II] 
122 ^tm data which solely traces this medium. In the 
absence of additional data, which would allow to seper- 
ate the contributions from the different phases of the 
ionized medium, we argue that dense H n regions are 
the primary source of [N n] emission. These emit 15% 
to 30% of the observed [Cn] emission while the re- 
mainder stems from PDRs. 

— The gas heating efficiency was calculated from the ra- 
tio of the two major gas coolants [O i] at 63 fim and 
[Cii]pdr at 158 fiui versus the total infrared intensity. 
The efficiency is low in the centers showing ratios of 
(1.6-2.1) 10 -3 , while the outer positions show higher 
ratios of (2.5-3.6) 10 -3 . The latter efficiencies lie at 
the high end of efficiencies observed by MKH01 in a 
sample of 60 unresolved normal galaxies. 

— We fitted the observed line intensity ratios 
[Cii] PD r/[Ci](1-0), [Ci](1-0)/CO(3-2), CO(3- 
2)/CO(l-0), [C ii] /CO (3-2), and, [O il(63) ACnW 
to the predictions of the PDR model of lKaufman et alJ 
(|1999^ . The best fits yield densities of 10 4 cm" 3 and 
FUV fields of ~ Go = 20 — 30 times the average 
interstellar, almost constant at the six positions 
studied here at 80" resolution. This finding may be 
a selection effect, since we selected the four outer 
positions for their high star forming activity. More 
observations of other less prominent positions are 
needed to study variation of physical parameters over 
the galaxy surfaces. 

— Filling factors vary significantly between the center po- 
sitions and the outer positions in both galaxies. The 
filling factor </>tjv derived from the FUV field calcu- 
lated from the observed total infrared intensity versus 
the fitted Go lies at 1 at the spiral arm positions, with 
a scatter of less than 20%. In contrast, the nuclei show 
ratios of ~ 4 — 5. The density contrasts within the 
emitting gas must be high, given the local densities of 
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10 4 cm -3 and the average densities, derived from CO 
column densities, of less than 1cm -3 , leading to very 
low volume filling factors. In accordance, the area fill- 
ing factors of the [C i] emission are less than 2% at all 
positions. 

— The reduced chi squared lie between 3 and 5 at three 
positions, including the two nuclei, while the fits are 
worse at the three other positions, which show \ 2 = 
13-15. It is shown that the fits can be significantly 
improved by assuming that [C i] emission is more ex- 
tended. 

While the present analysis led to the results listed 
above, it needs to be refined with new data and improved 
modelling. The assumption that all gas tracers have the 
same filling factors may not be justified. We find indica- 
tions that [Ci] is more extended than CO 2-1. Clearly, 
maps of [Ci] at 10" resolution filling the ISO/LWS 80" 
beam would be needed to verify this conclusion. A re- 
lated question is the nature of the interclump medium 
and whether it can be ignored when interpreting the 
submm/FIR emission. The fraction of [Cn] emission from 
the different phases of the ionized medium has been ad- 
dressed here by assuming that all emission stems from 
the dense ionized medium. While this is consistent with 
the PDR modelling, it probably is an oversimplification. 
Observations of the [Nil] line at 205 ^m would allow to 
refine this analysis. Such observations will become possi- 
ble with a new generation of submm telescopes and re- 
ceivers operating at high altitudes. In this regard, mod- 
elling neads to include the ionized medium. Improved 
PDR models need to consider a distribution of clouds 
within each beam, following mass and size distribution 
laws. 
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Fig. 4. Comparison of the observed line intensity ratios [C ii] PD r/[C i](l-0), [C i](l-0)/CO(3-2), 00 (3-21/00(1-0), 
[Cii1pdr/CO(3-2) and [Oi](63 itm)/[Cii] PD R, in M83 and M51 with PDR model calculations bv Kaufman et alJ 
(1999). The drawn contours correspond to the observed intensity ratios while the dashed contours are those for the 
20% uncertainty. The greyscale images show the reduced x 2 of the fit. The position of the minimum reduced chi 
squared is marked by a white cross (TableUJ) . The horizontal black line shows the FUV flux Go,obs calculated from the 
observed TIR C flux densities. 



